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THE CALCULATION OF DOWNWASH BEHIND SUPERSONIC WINGS 
WITH AN APPLICATION TO TRIANGUIAR PLAN FORMS 
By Max. A. Heaslet and Harvard Lomax 


SUMMARY 

A method is developed consistent with the assumptions of small 
perturbation theory, which provides a means of determining the down- 
wash "behind a wing in supersonic flow for a known load distribution. 
The analysis is based upon the use of supersonic doublets which are 
distributed over the plan form and wake of the wing in a manner 
determined from the wing loading. 

The first application of the method proves the equivalence in 
subsonic and supersonic flow of the dovnvash at infinity correspond- 
ing to a given load distribution. The principal application in thiB 
report is concerned with the downwash behind a triangular wing with 
leading edges swept back of the Mach cone from the vertex. A complete 
solution is given along the center line of the wake and an approxima- 
tion provided for points in the vicinity of this line. 


INTRODUCTION 

The linearization of the partial differential equation satisfied 
by the velocity potential for compressible flow yields, for subsonic 
flight speeds, an elliptic-type equation which is reducible by means 
of an elementary transformation to the basic equation in incompressible 
flow. As a consequence of this result, wing theory in the subsonic 
realm employs the same concepts and same types of analyses that 
belong to classical incompressible theory. At supersonic speeds, the 
differential equation for the velocity potential is hyperbolic in type 
and for wing theory is equivalent mathematically to the two-dimensional 
wave equation of physics. In spite of the different character of the 
basic differential equation in the two flight regimes, certain formal 
equivalencies can be set up which are intuitively useful in the solu- 
tion of specific problems. In particular, the velocity potentials 
of a three— d imens ional source and of a doublet each have analogous 
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forms in the two cases. The solution of different boundary—' value 
problems encountered in wing theory has been discussed in reference 
1 , and it has been shown how suitable distributions of sources and 
doublets may be used to determine the flow potential associated 
with a given lifting or nonlifting wing. 

The calculation of downwash behind a wing, for incompressible 
flow, relies almost exclusively on the use of Prandtl * s lifting— line 
theory which is, in turn, developed from the concept of a single 
horseshoe vortex. The conventional approach to the general down- 
wash problem is to determine, first, the induced field of the simple 
horseshoe vortex by means of the Biot-Savart law and, then, from a 
knowledge of the spanwise distribution of loading over the wing, to 
calculate finally the induced field produced by a vortex sheet 
composed of superimposed vortices of varying span. 

When downwash calculations are to be extended to the case of 
supersonic wings, it appears at first that the use of vortex Bheets 
is completely inadmissible since no equivalent to the Biot-Savart 
law exists. It is, in fact, true that the horseshoe vortex no 
longer plays the outstanding role it has at low speeds. However, 
when a more detailed investigation is made of the underlying analysis, 
it becomes apparent that vortex theory and the Biot-Savart law can 
be developed from the initial use of a constant distribution of 
doublets over a given surface, (l.g., see references 2 and 3 .) These 
doublets produce a discontinuity in the velocity potential at the 
surface, and, for incompressible theory, the curve which bounds the 
surface can be identified with a vortex curve possessing circulation. 
The proof of the Biot-Savart law and the introduction of vortex sheets 
are direct consequences of these basic ideas. 

Since, as was shown in reference 1, supersonic boundary— value 
problems involving sources, sinks, and doublets can be solved in a 
manner analogous to that used in low— speed theory, a method is 
therefore provided for an attack on the downwash problem for 
supersonic plan forms through the use of doublet distributions. 

The present report has two principal aims: First,' to outline 

in some detail the theoretical approach to the determination of the 
velocity potential of the flow field associated with a supersonic 
lifting surface and the subsequent calculation of the downwash; and, 
second, to apply the theory to the case of a triangular wing swept 
back of the Mach cone and to present the results of the complete 
calculations along the center line and in the plane of the wing. 
Downwash immediately back of the trailing edge and at an infinite 
distance behind the wing will also be derived and is in agreement 



NA.CA TN No. 1620 


3 


with the previously published results of P.A. Lagerstrom (reference 
4) . The final portion of the analysis will treat the variation of 
downwash near the center line of the triangular wing. 

In the theoretical portion of the report* the boundary— value 
problem will be introduced and the solutions* obtained from Green’s 
theorem* will be given for low— speed and supersonic flow. In the 
section of the report devoted to applications* the theory will be 
used first to evaluate the potential function at an infinite 
distance downstream from a lifting wing. Since the mathematical 
problems arising in the physically obvious case of the unswept wing 
of infinite span correspond closely to those for the more general 
case* the theory is next applied to this case. From this applica- 
tion, a general procedure is developed for treating wings with 
supersonic trailing edges. The final application of the report 
will be devoted to the triangular wing. In all of these applications* 
it will be seen that the analytic expressions which have been obtained 
in supersonic theory for the load distributions over certain plan 
forms afford a means whereby the chordwiBe distribution of pressure 
may be introduced into the analysis* and* therefore* such expedients 
as lifting— line theory are no longer so essential. 

No attempt will be made here to discuss the effect of airfoil 
thickness on the downwash distribution, although this effect is 
actually simpler to treat mathematically. It suffices to state 
that the entire theory is postulated on the assumptions of thin— 
airfoil or small— perturbation theory and that, consequently, thick- 
ness effects and lifting-plate solutions are additive . For the 
results that are given in the plane of the airfoil, the thickness 
effect* which Is necessarily symmetrical with respect to this plane, 
is zero. 


LIST OF IMPOBTANT SYMBOLS 


a Q velocity of sound in the free stream 

b wing span 


c 


o 


root chord of wing 




H 


oamplete elliptic integrals of the seoond kind with 
moduli k* k 0 , ki* ks * respectively 

gq,v 0 

s 0 P 
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x 0 -l 

Xq-1 


complete elliptic integrals of. the first- kind with 
moduli kj k 2j respectively 


free-etream Mach number 


static pressure 


ft) 


Pr-Pu 

dynamic pressure 


(|PoV) 


'( x-x x ) 2 + (y-yi ) 2 +( z-z ! ) 2 

f (x-x x ) 2 +(y-yi ) 2 +z 2 
k x-xi ) 2 -p 2 ( y-y i) 2 -p 2 ( z -z i.) 2 


'(x-xi) 2 -p E (y-yi) 2 -p 2 z 2 


u,t,w perturbation velocity components in the direction 

of the x—, y— , z— axe s , respectively 


x,y,z 


v-u-v-l 

free— stream velocity 

z— component of velocity induced by doublet distribution 
over plan form 

z— component of velocity induced by doublet distribution 
over wake 

Cartesian coordinates of an arbitrary point 


x 1j7ij z i Cartesian coordinates of source or doublet position 
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7o 

Zo 

a 

P 

6 

11 

So 

V 

Po 

4 

A4 S 


x/c 0 
P y/c 0 
Pz/c 0 

angle of attack 

Vm 2 -i 


semivertex angle of triangular wing 

PZi 


c o 0 o 


P tan 8 

Mach, angle ^arc sin ^ ^ 
density in free stream 
perturbation velocity potential 

$u~®Z 


sign denoting finite part of integral (equations 
(10) and (11) ) 


Subscripts 

u conditions on upper portion of surface 

l conditions on lover portion of surface 

L.E. conditions at leading edge 

T.E. conditions at trailing edge 

W refers to vake 

P refers to plan form 

s conditions on discontinuity surface (at zx=0) 
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Regions 

I, II, III integration regions on plan form (fig. 4) 
A,B regions in wake of triangular wing (fig. 5 ) 


THEORY 

Boundary Conditions 

Hie proposed problem is one of finding the downwash behind a 
flat plate which supports a loading consistent with its angle of 
attack and plan form. It will be assumed throughout the analysis 
that this load distribution is known. Such values were given for 
several plan forms in reference 5 and further results can be found 
in the literature on supersonic wings. 

The load distribution over the wing may be obtained from a knowl- 
edge of the differences in pressures acting on the lower and upper 
surfaces. Moreover , in thin-airfoil theory, where boundary conditions 
are given In the z = 0 plane (i.e., the plane of the wing), a simple 
relation exists between local— pressure coefficient and the streamwlse 
component of the perturbation velocity. Thus, assuming that the free— 
stream direction coincides with the positive x-axis (fig. 1), and 
denoting by u the x-component of the perturbation velocity, it 
follows that 

An p,— p 2 2Au 

— - hi. _ ( Uu ^,) a a) 

9 1 V 0 - v 0 


where the variables are defined in the table of the symbols. 
Furthermore, from the definition of the perturbation velocity poten- 
tial <j> -■ 


$ = / udx (2) 

where a is a point in a region at which the potential is zero. 
Combining equations (l) and (2), the jump in potential in the plane 
z=0 can be determined by integrating the Jump in the u induced 
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■velocity or, what amounts to the same thing* 
cient. Thus, 


px y p ' 2 - 

A * s= Vl.E. 


where the integration extends from the leading edge to the point x 
and A$s represents the jump in ® In the xy-plane . Since load 
coefficient Ap/q must be zero off the wing and since u is an odd 
function in z, the value of u must be zero for all points off 
the wing in the xy— plane . It follows that A$ a remains constant 
at a given span station for all values of x beyond the trailing— 
edge position. 

Figure 1 indicates an arbitrary lifting surface in the z=0 
plane together with the distribution of A<E S for given constant 
values of y and x. In both subsonic and supersonic theory, the 
wing together with the semi— inf inite strip extending downstream of 
the wing form a discontinuity surface for the velocity potential, 
while A5 > b is equal to 0 throughout the 'remaining portion of the 
xy— plane. These conditions, together with the fact that the - vertical 
induced velocity w Is a continuous function at z=0, are sufficient 
to determine <X> throughout space. The values of u, v, and w can then 
be found from the corresponding partial, derivatives of $ with respect 
to x, y, and z. The attention In the present report Is centered on 
w, the downwash f unction. 


the jump in load coeffi— 

(7) t3) 


Solution to Boundary— Value Problem 

In reference 1, the solutions for b oundary— value problems of 
the type under consideration were given for both incompressible and 
supersonic theory. The basic differential equations satisfied by 
the perturbation velocity potential are, for the two cases, respectively. 


d 2 I> 

dx 2 By 2 3z 2 


w 


0 
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and 


P 2 


d 2 $> 

8x 2 




15) 


Incompressible theory .— For boundary conditions prescribed in 
the z = 0 plane; the solution of equation (4) is 


<Sf (x^z) 



( 6 > 


where 


r = J { x-Xi ) 2 +(y— yi ) 2 +( z-z x ) 2 


and t is the area for which the integrand does not vanish. 

1 and l/8 

4irr e kn \3z! ±/ B 

x,y,z of a unit source and doublet situated at the point 


terms 


- 

: i ry £ 


The 


are equal to the velocity potential 


-s 

at x,y,z of 

z ljyij 0. The remaining terns in the integrand, which determine the 
distribution of source and doublet strengths, must be found from 
known boundary conditions. If a lifting surface fixes the boundary 
conditions, induced vertical velocities on the upper and lower 
faces of the surface are equal so that 


5ft i 

Szi dzi 


and 

* (w) ' kf r ST £). 31145,1 l7> 

Equations (6) and ( 7 ) are well known in potential theory (refer- 
ence 3 j P. 60), but the derivation usually employs the assumption 
that the value of ® is zero, at all points infinitely distant from 
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the wing. This assumption cannot , of course, be made in aerodynamic 
applications where the discontinuity surface t extends to x = <*>, 
as in the case of a lifting wing or lifting line with trailing 
vortices. These latter problems, with which this report is directly 
concerned, are of such a nature, however, that the induced effects 
at an infinite distance are confined to the plane x = 00 . An investi- 
gation of the derivation of equation (6) reveals that the conditions 
imposed on $, in general, can be relaxed sufficiently to permit a 
discontinuity in a strip of finite width along the entire extent of 
the x-axis. The mathematical details of the derivation will not be 
given here but a statement of the restrictions on 4 at infinity is 

worthwhile. Thus* denoting hy — the directional derivative of 0 

dn 

taken normal to a prescribed surface, the following conditions 
apply: 

d5> 

1. The functions <!> and - — are zero at all points 

on 

having radius vectors which make finite (nonzero) angles- 
with the positive x-axis, the points lying on a spherical 
surface of infinite radius with center at the wing. (This 
preserves the usual potential theory assumptions except 
over the portion of the spherical surface -which forms the 
plane x = oo. ) 


2. The values 0 and r— are bounded at all points 

ox 

infinitely distant from the lifting surface and at a non- 
infinite distanoe from the positive x-axis. (This condition 

places restrictions on the values of $ and — in the 
plane x = <=o.J 


Conditions (l) and (2) are satisfied for a lifting surface of 
finite span, and equation ( 7 ) is consequently applicable' directly 
to the determination of the velocity potential. As an application 
of the equation, suppose a sheet of horseshoe vortices is situated 
as in figure 2 with bound vortices placed on the y-axis, trailing 
vortices extending parallel to the positive x-axis, and has a span- 
wise distribution of circulation A$ symmetrical to the xz-plane 

"b "b 

and defined for — < y < ^-. Then the velocity potential corre- 

sponding to this vortex sheet is given by the expression 




dxi 


' 0 [ (x-xi ) 2 +(y-yi) 2 +z 2 3 3 2 


( 8 ) 
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When A$ s = constant, a single horseshoe vortex results. 

Supersonic theory .— For supersonic boundary— value problems 
associated with plan forms as indicated in figure 1(a), where the 
known conditions are given in the z = 0 plane, the general solution 
of equation (5) is given in reference (1) in the form 


® (x,y,z) 



where 


r 0 = v /(x-x 1 ) 2 - 13 2 (y-y 1 ) 2 -p 2 (z-z 1 ) k 


and the subscript s on the parenthes-es indicates that the function 
is to be evaluated at z x = 0. The region t is that portion of the 
Xxyi— plane bounded by the leading edge of the wing, the lines parallel 
to the x— axis and stemming from the lateral tips of the wing, and the 
trace in the z x = 0 plane of the Mach forecone with vertex at the 

point x,y,z. The sign [ is to be read "finite part of" 

and was introduced by Hadamard (reference 6) as a manipulative 
technique with the property that 



A(x)dx _ rx o A(x) — A(xq ) 
(x 0 -x) 3/2 J a (x 0 -x) 3/2 


dx — 


2A(*o) 

(xo-a) 17 ^ 


( 10 ) 


For purposes of calculations, this was modified in reference 1 to 



A(x)dx 
(x 0 -x) 3 / 2 



A(x)dx 
(Xo-x) 3 / 2 


= - F (a).- C 


( 11 ) 


the asterisk indicating that no upper limit is to be substituted 
into the indefinite integral, the latter being determined as 


F (x) + C 
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where 


C = lim 
x x 0 


- , w ‘ 


Equation ( 9 ) is the direct analogue of equation (6). The terms 

and ^ ~ are equal to the velocity potential at 

x,y,z of a unit supersonic source and doublet situated at the point 
while the remaining terms in the integrand determine the 
distribution of source and doublet strength and are determined by the 
known boundary conditions. 

When the potential function associated with a lifting surface is 
to be evaluated. 


1_ 
2 it 


c/ s 



Szi Szi 


and equation ( 9 ) reduces to the form 


® 


IJ (sr £> 4X1411 


( 12 ) 


In application, the region of integration in equations ( 7 ) and (12) 
can be divided into areas occupied, respectively, by the plan form 
and the wake region. Thus, for equation (12), 


® (x,y,z) 


-zp g r r A3> s dx^dyj, 

2it J J [(z-x 1 ) 2 -p 2 (y-yi) 2 -P 2 z 2 ] 3 / 2 

plan form 


zP 2 

r r 

A® s dXidyx 

2 n 

^ wake 

[(x-xi ) 2 -P 2 (y-yi ) 2 -P 2 z 2 ] 3 / 2 


( 13 ) 
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Equation (13) presents a formal solution for the calculation 
of velocity potential and, subsequently, downwash for a given surface 
in terms of A4 B . Since A$ e was related directly to load distri- 
bution in equation (3), it is apparent that the various known 
solutions to lifting— surface problems are directly applicable. The 
fact that supersonic theory permits the determination of load 
distribution in closed analytic form for many simple plan forms 
provides a distinct advantage that is lacking in subsonic theory 
wherein virtually all known results are available only in numerical 
form. Thus, theoretioal analysis of problems involving supersonic 
flight speeds can be carried further before recourse to numerical 
methods is necessary. 


APPLICATIONS 

Value of Potential Function at x = » 

It is possible to show, from equations (7) and (12), that the 
potential functions corresponding to a wing with fixed load distri- 
bution are idential at- x = » for incompressible and supersonic 
flow. Assuming At> a known, the values of $ (x,y,z) for the two 
cases are given, respectively, by the equations 


® (^y,z) 


z nr dfl> B dxidyi 

J J E (x-xi ) 2 +(y-yi ) 2 +z?l 

plan form 


*hff 

wake 


A$s dXjdyi 

[(x-xi ) 2 +(y-yi) a +z 2 W s 
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®(x,y,z) = 


plan form 


M b dxxd y x 

[ (x-xx) a -g 2 (y-y 1 ) 2 _p 2 z2] 3/2 


-S-UT 


A$ g dxxdyx 

[(x-Xx^^y-yx) 2 -^ 2 ] 3 / 2 


Since, however, A$ s is finite, it follows Imm ediately that for 
fixed values of y and z the integrals over the plan form in "both 
equations approach zero as x increases indefinitely. Thus, 
denoting "by x^.-g. the value of xx at the trailing edge of the 
wing, the potential functions at x = 00 are given by the expressions 


vSt. hf***^ 1 *-*-*^/ 


XT.E. 


[ (x-Xx) 2 +(y-yx) 2 +z 2 ] 3/2 


«(“* 7,z) = 11m — — f . ADgC^ E ,y )dy f 
x-> <» 2 jt >jL1l J 

2 


^T.E. 


[ (x-xx ) 2 HB 2 (y-yx ) 2 -P 2 z 2 ] 3 / 2 


These relations can be integrated once to give for the subsonic case 


O(co,y,z) = lim ^ r&bJxrn -m ,y )dy, "I 

X co J §. * l [(y-yi ) 2 +z 2 l / (X-Xx ) 2 +(y-yx ) 2 +z 2 J 


1_ Pz A^s 

2rtJb (y-7i) 2 +z 2 


(14a) 



and for the supers cnio case 


£ 


h * 

$U,y,z) = lira. — P A® a ,y, ) dy. X ... , 

I- >* 2s ^ ' l Ky-yi) +z a 3 i/(x~x 1 ) E -p 2 (y-yi)^-p 2 z 2 Ji T< j > 


fc. 

!_ f 2 

at J.t (y-yi) a +z 2 


(ito) 


From these equations, it follows that the side wash and downwash at x = ■» are invariant 
with Mach nuinbe^ provided the load distribution is fixed. In fact, their talues depend solely 
on the spanwise load distribution, since the terms corresponding to the ohordwise distri- 
bution disappeared in the analysis. This has been pointed out elsewhere in the literature. 

It should be stressed, however, that the result which has been obtained here states that 
equal span load distributions in the two cases yield equal values of the potential function 
at X ^ to. This does not imply that a wing at low and supersonic speeds maintains the same 
potential function at Infinity. When the wing is kept fixed, the distribution of w on 
the wing is fixed, but the load distribution la a function of speed. 


Downwash On and Off the Wing 

As a further application, the unswept wing of infinite span will be treated for superb- 
sonic speeds. In this case, as is well known, the induced velocities are zero at all points 
downstream of the upper and lower Mach waves stemming frcrn the trailing edge. An abrupt 
Jump in vertical velocity therefore occurs at the trailing edge of the wing. Consideration 
of this jump for the unswept wing furnishes considerable insight into the nature of the 
mathematical difficulties inherent in the calculation of downwash on and off wings of 
arbitrary plan form. The calculation for the particular case will therefore be followed by 
a more general discussion which will be of value in connection with the later treatment of 
the triangular wing. 
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Uhswept wing of Infinite span .- The pressure distribution for the wing of infinite 

aspect ratio is constant. For this so-called Ackeret-type loading , ^ is eq ual to — 
so that, when the leading edge lies along the y-axis , ^ ^ 


A4 C 


Vo r x ^ 

2 1 q 


dx.^2 

5 P 


*1 


( 15 a) 


where a is angle of attack. In the wake 

,°o 


A<$„ 


_lo r 
2 J q 


2aV 0 
dx = — 

q P 


(15b) 


The downwash, or vertical velocity, will first be found when the point is between the Mach 
waves from the leading and trailing edges and then when the point is downstream of the 
trailing-edge wave. 

Since the wing is of infinite width and eafperiences no variation with y, it is 
possible to consider the problem at y = 0. Thus, from equation (13), for the case when 
the point under consideration is between the Mach waves from the leading and trailing 
edges with its forecone cutting the wing as shown in figure 3(a), 


. -zf^Voa 

= lxm 

e o tcP - 0 


/ 


x— j3 ./e 2 +z 2 


xi dxx 


L 


^V(x-xi) 2 -pV 


iyi 


ft x-xx ) 2 H3 2 yx 2 -p 2 z 2 ] 3/2 


H 

VJ1 
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where the symmetry of the problem with respect to the yx-axia has been used. Computing the h 
finite part of the integral ^ 


= lim 
e -^o 


g2py 0 a g r x ~Q V^ 2 

It Jo 


*1 dx* 


[(x-Xx ) 2 -? 2 ^ 2 ] y ( x-xx ) s -& 2 e 2 -£ 2 


(16) 


and 


8$ 5 255p?oae ^ 2+z2 

w = — = — lim / 

Sz bz e - >' 0 it J n 


Xi clxi 


Kx - xx ) 2 -^^ 2 ] yu - xi ) 2 -^ 2 -^ 2 


( 17 ) 


Equations (l 6 ) and (17) can be evaluated directly to give the results 


V n <x z 

4> * (x±pz) 

0 z 


(18) 


v = — 


V 
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For a point behind the trailing-edge wave (fig. 3(b)), the two quantities can be 
determined in a similar manner. Thus, 


fl = lim 

€ 0 


-gzpYoae P G o Xi dxx 

it J Q [(x-xx) 2 -p 2 z 2 ] / (x—xx ) 2 -p 2 z 2 -p^£ 2 


2zpYoac #-^) £ -P 2 z 2 
+ / p dy x 


/ 


x-p /yi 2 +z 2 


dxi 


[(x-Xx^yi 2 -^ 2 ] 3 / 2 


(19) 


and v iB given by the partial derivative of $ with respect to z. The term containing 
the single integral ia zero, since the integral itself is bounded for all values of e, 
while the term containing the double integral is readily calculable. Thus, the values of 
the velocity potential and the dovnwash for a point behind the tralllng-edge wave are given 
by the relations 


<P = - 


YpOCo 

p 



w = 0 


( 20 ) 


These results are the familiar equations associated with two-dimensional supersonic flat- 
plate theory. 


f— 1 

-4 
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The point of principal interest in this development is the 
jump in the induced vertical velocity w in passing from a point 
just ahead of the trailing-edge wave ter a point immediately behind 
the wave. A study of equations (17) and- (19) shows that this Jump 
is the result of the discontinuity in the contribution to the down- 
wash of the term containing the single integral. Ahead of the 
trailing-edge wave, this term yielded the result that 


w = — Y <x 
o 

whereas behind the wave, the contribution of the term to w was 
zero. 


The method of attack used in the study of the unewept wing can 
be generalized to apply to arbitrary plan forms. A discussion of 
this case follows. 

Arbitrary pTa-n forms .— As will be shown later, for any plan 
form with supersonic trailing edge, the jump in the value of w 
in the plane of the wing at the trailing edge can be calculated 
directly by means of simple momentum methods. At this point, how- 
ever, it is of more interest to consider in a genral manner the 
nature of the integrations involved when the point x,y,z is either 
ahead of or behind the trailing-edge wave. Figures 4(a) and 4(b) 
show a plan form with a straight trailing edge with areas of integration 
indicated for the point P in each of the two positions. (The 
straight trailing edge is not a necessary restriction and is only 
introduced for convenience of notation.) The regions of integration 
are divided under the assumption that the first integration on the 
plan form in equation ( 13 ) will be made with respect to y x . When 
the point P is ahead of the trailing-edge wave, therefore, the 
contribution of the wake is zero and the integration over the plan 
form is made to conform with regions I and II. When the point is 
behind the trailing-edge wavey three integrals are evaluated 
corresponding to regions I, II, and HI. In the case of the infinite 
aspect ratio, unswept wing region I was, of course, nonexistent and, 
in general, no essential difficulty in regard to the limits of 
integration Is introduced by this region regardless of where P Is 
situated. In region II, however, the problem must be treated in 
more detail. 

Consider first the case when P is ahead of the wave and denote 
by $Ha the contribution of one side of region II to the total 
potential. Then 
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Similarly, when P is 'behind, the wave and the same subscript 
notation is used to refer to one side of region II, the value of 
4 Ila is 


®Ha 


= lim 

e -3»- o 


-ft 22 r 00 ^ r Yl ASi UijyiJdyi 

2* Jx 1 ^ y+e [(x-xi) 2 -p 2 (y-y 1 ) 2 -p 2 z 2 ] 3 / 2 


(21b) 


where Tx is as defined above. 

The contributions of the other side of region II to the potential 
will not be considered separately as the behavior is identical. When 
P lies ahead of the wave, e appears in the limits for integration 
with respect to both j x and s*. This corresponds to the situation 
in equations (l6) and ( 17 ) and, as for that problem, the limiting 
process is carried out after the integration is completed. When P 
is behind the wave, it is not necessary to defer the limiting process, 
since 



t* 
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But since 


ggi (xi,yi) ' 

[(x-xO^^y-yx) 2 -^ 2 } 3 /* 


< M 


for Co = Xi = Xo (i.e., is Bounded for all values of xx in the 

interval of the first integration); and further, since 


r? +£ 

lim M / dyi = lim Me = 0 
e o J y s 0 


therefore, for P situated Behind the trailing-edge wave, the 
contribution of region Ila is given By 



(x 3 .,yi)dy 1 

[(x-Xx^-p^y-yx )^^ 2 ] 3 / 2 


( 22 ) 


Equation (22) will Be applied directly in the determination of 
downwash behind the triangular wing. The significance of the result 
is that, when the point' P at x,y, z is behind the Mach wave from 
a supersonic trailing edge, the limiting process associated with 
region II need not Be considered. When P is ahead of the Mach wave, 
the term e must Be retained in the analysis and the limiting process 
used. As was previously noted, the general analysis developed in 
this report places no restriction on the orientation of the trailing 
edge; however, it should Be pointed out that region II exists only 
for the case in which the trailing edge is supersonic. Therefore, 
the Jump in downwash, obtained from the integration over region II, 
is associated only with supersonic trailing edges; whereas Both the 
downwash and loading are continuous across a subsonic trailing edge. 


Triangular Wing 

Consider a triangular wing (fig. 5) with leading edges swept 
Back of the Mach cone from the vertex. The loading over the wing 
is known to Be (references 7 and 5 ) 
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Ap US 0 2 cc^i 

q E 0 P t/eo 2 ^! 2 -^ 2 ^! 2 


(23) 


where E 0 is the complete elliptic integral of the second kind with 

modulus k 0 = f/l—S Q 2 and 6 a = J3 tan S, 6 heing the semi vertex 
angle of the triangle. From equation ( 3 ) 


A3» s = H V’So 2 *! 2 — P 2 y:i _ 2 


(24) 


where 


E 0 P 


(25) 


Sufficient information is now at hand to permit the use of equa- 
tion (13). Setting, for convenience, ®= <&p+ 3> w the velocity 
potential at x,y,z is given hy the sum of the two expressions 


Gp = — 


zH[3‘ 

2it 


IL 


■/go 2 ^! 2 -^^! 2 dx x dy x 
plan form [(x-Xi) 2 -p 2 (y-yi) 2 -p 2 z 2 ] 3 / 2 


( 26 ) 


®W = “ 


zHP‘ 

2k 


//„ 


wake 


//e 0 2 c 0 2 -p 2 yi 2 dx x dy x 
[ ( x— xj. ) 2 — {3 2 (y— y 1 ) 2 — p 2 z 2 ] 3 / 2 


(27) 


Equation (26) represents the contribution to the velocity potential 
furnished by the doublets distributed over the plan form, while 
equation ( 27 ) represents the contribution furnished by the doublets 
within the wake . The latter equation is the mathematical equivalent 
in supersonic flow of the subsonic velocity potential of a sheet of 
horseshoe vortices corresponding to an elliptic span load distri- 
bution. Equations (l4a) and (l4b) showed that the expression for 
at x - co is identical to the velocity potential of the 
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subsonic vortex sheet. However, in the vicinity of the x = c 0 
line, the behavior is entirely different. 

In the present report, equations (26) and ( 27 ) will be applied 
only to the ds'tezminati on of downwash along the center line in the 
wake of the airfoil. Further extensions are limited only by 
difficulty of integration. It will appear in the development that 
recourse to numerical methods is probably necessary in general. 

Downwash induced by doublets in the wake .— Setting y = 0 in 
equation ( 27 ) and integrating with respect to Xi, it’ follows that 


zH(x-cq) V / e 0 a Co g -P g yi g dy x 

2lt -Yi (yi 2 +z a ) V( x_c o ) g -P g yi a -£ g z g 


( 28 ) 


where the limits on the integral are not yet specified, since they 
differ in the regions A and B indicated in figure 5- In either 
case, however, the limits are seen to be the roots of one of the two 
radicals in the integrand. 

In order to derive an expression for downwash in the plane of 
the airfoil, it is convenient to express equation (28) in a different 
form. Integrating by parts 


H 


ye^co 2 -^! 2 

2 * 


arc tan 


yiU-Cp) 

z y(x-c 0 ) 2 -£ 2 yi 2 -P 2 z a 



+ 


B 2 

2rt 



yi 


ye 0 2 c 0 2 — 


arc tan 


yi( x ~ c 0 ) 

z V(x-c o ) 2 -0 2 yi 2 -^ 2 z 2 


<3yi 


This form of the equation is easier to treat when lim - — is to 

z -> o oz 

be considered. In both regions A and B, it can be shown that w^, 
the contribution to the downwash made by the doublets in the wake, 
is given by the expression 
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v ¥ = 


HP 2 

2«(x— Co) 


■ Yl / U-ool^Vi* 

eo^o 2 ^ 2 ^! 2 


£7 


d-7i 


(29) 


In region B, . Yx = 


Qq c o 

T 


Using the fact that the integrand is 


an even function and introducing the notation 



(30) 


Setting kx = — — 

Xq-1 

expression for w^- 


and noting that kx < 1 in the region, the 
in region B is 


wvr = - 52- lx 


(3D 


where Ex is the complete elliptic integral of the second kind and 
the subscript denotes that the modulus is kx. 

In region A 


*1 


x-Qq 

P 


By means of the transformations 


x 



T) = SZi_ 

x-c D 


X 



2k 


M.QA T TS I fe* 1629 


the integral is transformed to the form 
v W 1t Pi 

HP' " ” Jn 


1 -n 4 




-&T) 


So— 1 1 

Setting k 2 = = =■, noting that kg < 1 in region A, and 

Q 0 kj. 

introducing Jacobian elliptic functions in the tra nsf orma t i on 
i) = m u, ve have 


Wu-Jt p Ke (l^k2 2 )Ee-$2 

= — k 2 / cn 2 udu = - _ 

Jo is 


nr 

HP 


(32) 


where K 2 and E 2 are complete elliptic integrals of the first and 
second kind with modulus k 2 . 

The values of the downwash given by equations (31) and (32) can 
be presented in terms of w 0 , the downwash cm the wing, where from 

/ \ ^-^Wq 

equation ( 25 ) and from the boundary condition, c t = =-=■, 

* n 


v 0 it 

w 


’ - 2 S ° 


Hus, 


w. 


_W 

Wo 


= 2 


Eg— ( tek 2 . 2 )Kg 
kE Q k 2 


(33) 


( region A) 


*W_ 

w 0 


2Ei 

w W 


( region 5) 


Figure 6 shows this variation in •= — plotted as a function of x Q 

V ° 

for values of 0 o equal to 0 . 2 , 0,4, 0 , 6 , and 0 , 8 . 

Downwash induced by doublets on the plan form .-r- Along the center 
line , the contribution to the potential fxmcticm duo to the distribution 
of doublets on the plan form is given by equation ( 26 ), where the value 
of y is zero, dust as in the case of the wake distribution, it is 
seen that the nature of the solution changes in passing from region A 
to region B 90 that tfe© derivations for t&©£© regions will b© given 
independently « 
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In region B, no ■ singularities occur in equation ( 2 6) so the 
finite part sign can he disregarded. Since , moreover, the integrand 
is an even function of y 1} the expression for velocity potential 
becomes 


9 n = - 


zH( 3 2 


gogi 

3 


r d^ rn dyi (31,.) 

Jo Jo [rx-Xi)2_p2 7l 2-p2 Z 2 ]3/2 71 


After changing variables by the transformations 


S = 


fiyi 

e 0 xi 


9q^i 

J (x-x x ) 2 -p 2 z 2 


equation (3^) becomes 


$p 





di 


Substituting |=sn u and noting that sn K=l, sn 0=0 where K is 
the complete elliptic integral with modulus k, we have 


<X>p = — 


zHp J C ° 


k 3 

e o X! 



cd 2 udu = 



(K— EOdxi! 

J (x-Xi) 2 -p 2 z 2 


The d<^mwash wp in the 

lim — r— . This leads to 
oz 


z=0 plane can now be found by considering 
the result 



(E-E) 

X— Xi 


dxi 


where the modulus reduces to 
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k 


X— Xi 


Rewriting the equation so that- k is the variable of integration 
and ki = 

Xq-1 


Wprt 

w 



(35) 


As is to be expected, the value of Wp is seen to approach 
zero as Xq becomes infinitely large. The upstream boundary of 
region B lies at Xo=l+£>o and, as will be seen, Wp is continuous 
at this point. 

Consider now region A. In this case, the traces of the Jfech 
forecone of the point (x,o,z) cut across the plan form. For ease 
of calculation, the area within the wing leading and trailing edges 
and the traces of the forecone will be divided into two parts. Thus, 
making use again of the symmetry with respect to the xi— axis, 
and also the results discussed in connection with equation (22), it 
follows from equation (26) that 



The reduction of the first double integral in equation (36) can be 
carried out by methods exactly equivalent to those used in reducing 
equation (3*0 to (35) ♦ As a result of this calculation, the contribu- 
tion made to Wp in the plane z=0 by the first double integral is 
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cUc 


(37) 


In tiie second double integral of equation ( 36 ) the following 
substitution is made 

I = Pyi _ //(x-x i) g -p a z a 

^/(x-xx) I 2 -^ a z 2 0 o zx 

Then., if 


I = 




j/U-xJ^z 2 /0 o 2 Xl 2 -p s yi 2 


[ (x— x 1 ) 2 -p 2 y 1 2 -p 2 z 2 ] 3 / 2 


cLyi 


I can he written 


I - 


1 T 1 
pcTPv (1-1 


x Vl-k 2 S 2 


peoXxk^J (1-r) 

o 


27372 


dl 


The integration can heat he completed here through the use of equa- 
tion (l4) where 


A(i) 


(l +|) 3 / 2 


and the auxiliary transformation g = sn u is used. As a result of 
this calculation 


I h— (K-B) 

P0 o xxk 2 

so that the second double integral of equation ( 36 ) is equal to 

zHp P c o E-E 

/ dxx 

« e 0 xxk 2 
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Using methods similar to those employed In the derivation of equa- 
tion ( 35 ) , the contribution to Wp is found to he 


II 

it 



K-E 

^(l+So*) 


dk 


(38) 


where ka =• . Equations (37) and ( 38 ) may now he combined to 

e o 

give the total downwash in region A induced by the doublets on the 
plan form. Thus, 


Wplt 

HP 



g-E 

^(l+So*) 



K-E 

k+Po 


dk 


(39) 


The integration of- equation (39) requires numerical methods. A 
slight simplification can, however, be introduced by using the known 
value of the downwash w 3 just behind the trailing edge. As will 
be shown in the next section. 


W*~i ***«' 


Equation ( 39 ) must, of course, yield this result for Xq on the 
trailing edge, that is, for Xq= 1 and k 2 = 0 . Using this relation 
together with equation ( 39 ) it follows that 


Wplt 

w 



K-E 

k*(ue 0 k) 


dk 


(to) 


Values given by equations ( 35 ) and (40) are consistent at the 
point ki=k 2 =l, that is, at the point where the Mach cones from the 
two trailing-edge tips intersect on the x— axis . By means of these 
equations, the contribution to the downwash of the doublets distri- 
buted over the plan form is determined. 
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In summary, these relations are 


= E 0 -9 0 _ _2_ p ** K-E 

w Q E 0 rtE 0 J 0 k 2 (l+e 0 k) 

2 , _L. r * 1 “ at 
w o *E 0 J 0 k_H? o 


(region A) 


(region B) 


w p 

Figure 7 shows the results of the Integrations, the function — 

w o 

Being plotted as a function of Xq for various values of d Q . 

Conditions at the trailing edge .— The value of the vertical 
induced velocity immediately ahead of and Behind the trailing-edge 
wave must, of course. Be determinaBle directly from equations (21) 
and (22), respectively. By setting x=c 0 ± z£. If, however, the 
discussion Is restricted to the z=0 plane, a much simpler method 
exists for finding w at these points. The approach taken here 
follows essentially that given By Lager strom in reference 4. 

Let conditions just ahead of the trailing-edge wave Be denoted 
By the suBscript 2 and conditions just Behind the wave By the 
suBscript 3. Figure 8 shows a section of a given wing in the plane 
y=constant . The Mach waves at the leading and trailing edge make the 

angle p=arc sin with the z =0 plane, and the wing is presumed 

to Be at angle of attack a. Assuming the trailing edge to Be 
normal to the free— stream direction, the variation in the x— component 
in velocity when passing through the trailing-edge wave can Be treated 
as a two-dimensional proBlem with the condition imposed that U 3=0 
In the z=0 plane. It Is known that continuity of flow together with 
Balance of tangential momentum across the wave lead to the result 
that the component of velocity tangential to the wave is continuous. 
The tangential components of velocity Im mediately ahead of and Behind 
the wave are given, respectively. By the expressions 


(T ^) 2 = (V 0 +U 2 ) cos + w 2 sin p 
( y t )3 = y o P + V 3 sin p 
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Equating these relations, it follows that 


w s = w 2 + u 2 cot p. 

< 41 > 

For the two plan f 037 ms which, have heen considered, equation (4l) 
gives the following ire suits 

(a) Unswept wing of infinite span 

I 21 this case 


(7) 

H 2 

and 

w 3 = - V 0 a + -j- — = 0 (42) 

Equation (20) showed that this result actually applies to all points 
behind the Mach wave . 

(b) Triangular wing 

From equation ( 23 ) 


4a, 

P 

oP 4ct, 


4e 0 g aCo 


\ q4 EoP //0o 2 co^-P 2 yi 2 


from which it follows that 


w 3 = 


v o ( 1 “ 


0 o 2 c o 


E 0 '/8o Z °o Z ~?> Z 7i Z 


) 


( 43 ) 
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Setting yi=0 and introducing H from equation (25) 

t£ = -f (io-eo) m 

which, is the egression used in equation (40). These results are 
equivalent to those given by Lagerstrom in reference 4. 

Values of downwash near center line of wake .— The values of 
downwash, which were obtained on the center line of the wake, were 
worked out in exact detail and subjected to no restrictions other 
than those originally imposed by the use of the linearized equations 
of flow. It is possible, moreover, to get an indication of the 
variation in the downwash function in a portion of the wake for 
points near the center line through the use of a generalized Taylor’s 
expansion in the vicinity of the line y=z=0. The next higher terms 
in the expansion can be found without too much difficulty for the 
region bounded as follows 

(a) - | b <y <| b 


(b) Both y and z lying within the Mach cones from the 
trail ing-edge tips 

The problem resolves itself into one of finding the first 
nonvanishing coefficients of y Q and z 0 in the series 

~ = A 0 +A 1 z 0 +B 1 y 0 -hA 2 z 0 2 +C 2 z 0 y 0 +B 2 y 0 2 +.. < . (h-3) 


The value of Aq has, of course, already been computed and is 
known from equations ( 31 ) and (35) "to be 



(k6) 


So . 

Xq-1* 


where kx = 
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The coefficient Bx in the expansion is known to he equal to 


(© 


The evaluation of this term follows a devlopment s imi lar 


y o = 0 


to that used for equations ( 29 ) and (34) except that y is retained 
in the analysis. The derivation yields the relation 


wit 1 P 9 o ^ /i 

HP ' 2(xq— l) Jq o ^ * 


(x 0 -i) 2 -(y 0 -n ) 2 

© n 2 -n 2 


dri 


- i r\t. r te ° rSlE 

2 J J ffv 


U 0 O c (x °“£ ) 2 -(y 0 ~ T i) 2 ) 3/2 


df] 


Carrying out the differentiation, with proper regard for the singu- 
larity in the first integral, it follows that Bx = 0. Similarly, 


it can he shown that C 2 = 0, while the coefficient B 2 
is given hy the expression 


= Ifd 2 * > 

2 V&y 2 / 


y Q =o 


b 2 = 


TtE 0 e 0 ‘ 


2— kx 

2Kx Ex 

1-kx 2 


(47) 


”®o (k+So) [ E 


1+ * 2 _r 1 


(l-k E ) : 


1 ~k Z . 


dk 


where again kx = 


x 0 -l’ 


In order to calculate the variation with z, it is necessary 
to evaluate Ax 


- f-^-) 


0 ' Z 0 =° 
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where 


— = P f 
HP J n 


6 0 C 0 

“T” 


7 1 


arc tan 


yi(x-co) 




<3yi 


- zS / tai / 


gpXl 

p f/o7^7-7j7 


[ (x-X!) 2 -P 2 7i 2 -P 2 z 2 ] 3/2 


d-yi 


The double integral contributes nothing to the coefficient, and the 
remaining portion of the expression can be evaluated without 
integrating by differentiating twice and using Cauchy’s integral 
theorem 


Thus , 



f (yi)dyi 
(yi-iz) 2 


2«i 



Ai = - 


1 z 
EqSq I z 1 


(48) 


The coefficient A 2 will not be evaluated, Bince the first 
higher order term in z has been found. Thus, to the first order 
in y Q and z Q , the downwash function — is 


w 2_ 

w o ^o 





Zp 

= 0*0 



(49) 


DISCUSS ION 

The variation of — , the downwash due to the doublets in the 
w o 

wake in terms of the downwash on the triangular wing, was given in 
figure 6 as a function of x 0 for various values of 6 0 . This 
distribution of downwash along the center line of the wake can be 
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presented as a single curve 
A graph of this function is 


if — E n is expressed in terms of 

v 0 

shown in figure 9 . 


x 0 -l 
So ' 


Figure 10 presents the resultant downwash on the center line 
induced by the doublets distributed over the wake and over the plan 
form. This is the complete, value of downwash on the XQ-axis and 
shows that, following the discontinuity in downwash at the trailing 

edge, the magnitude of — builds up to its asymptotic value and 

v o 

achieves an almost constant value within approximately one chord 
length aft of the trailing edge for values of 0 O less than 0.8. 

The table given on the figure relates 9 0 to free— stream Mach number 
M and to 0 for the particular case of a triangular wing with lead- 
ing edges swept back 45°. 


It has been shown that at x = » the downwash has the same form 
for both incompressible and supersonic flow, provided the span load 
distributions are equal. An exact solution for the downwash on the 
center line of the triangular wing has been developed and extended 
to include a region near the center line of the wake by means of an 
approximation (equation (45) ) . A measure of the exactness of the 
approximation may be obtained by comparing these results with the 
exact results for the incompressible case for an elliptic span load 
distribution, as obtained from reference 8 , page 151. Figure 11 
shows the exact values of JL for incompressible theory for 

^ _ w o 


o< IS. 

9 o 


<0.3 and 0 < — < 0.3. 

9 o 

mations in supersonic . theory to 


w 

w. 


Also included are the linear approxi— 
for the same range of the 


E r 


geometric variables. The span load distributions for both wings are 
elliptical, but, in order that the same value of lift is maintained 
for the low— and high-speed cases, the angle of attack must be 
modified. This accounts for the use of the factor E Q in one apt 
of ordinates. It is apparent from the comparison that, at a large 
distance behind the wing, the approximation of equation (49) yields 
results within 4 percent of the exact value for the region considered. 


In application, it is desirable to know not only the order of 
the variables retained in the analysis, but also the magnitudes of 
the neglected terms. In this respect, equation (49) does not supply 
as. .much Information as might be desired, since no estimation Is 
furnished of further coefficients in the series. It appears, however, 
from the nature of the agreement obtained in figure 11 for conditions 
at a large distance from the wing together with observations as to 
the asymptotic behavior apparent in figure 10 , that equation (49) 
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furnishes a reasonable estimate of downwash in region B, provided 
that the point of calculation is not in the immediate vicinity of 
the tip Mach cones. For values in region A, no conclusions can he 
drawn from this investigation about the variation of downwash with 
respect to z. The spanwise variation indicated by equation (k 9 ) 
agrees with the value given by equation (V3) at the trailing edge 3 
that is, to the first order in y, the downwash is a constant. 

Thus, in the xy— plane, equation ( 49 ) can be expected to give a 
reasonably close estimate to the downwash near the center line and 
away from the tip Mach cones for all points behind the trailing edge . 
For variation in the z-direction, however, the validity of equation (^ 9 ) 
is restricted to. points in region B. 

The behavior of the downwash in the vicinity of the trailing 
edge indicates that in this region an important difference exists 
between lifting— line theory in supersonic and in subsonic flow. 

Figure 6 shows the variation of downwash along the center line of 
a doublet sheet (lifting— line theory in incompressible flow) with an 
elliptic span load distribution. The error introduced by not using 
the chordwise distribution of load is given in figure 7 and is 
sizeable for about one chord length back of the wing. 

The methods of analysis presented here were shown in the first 
application to lead immediately to the value of the potential 
function at x = 00, and thus provide a ready method for the deter- 
mination of vortex drag of a supersonic wing. 

A more detailed study of downwash behind the wing will necessarily 
involve considerable labor. The methods given in the report are, 
however, general and are directly applicable. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif. 
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( a ) Plan form 




(b) Sections showing distribution of A$ s 



Figure /. - Sketch showing arbitrary lifting surface together 
with distribution of A$ s , the jump in perturbation velocity 
potential in the plane of the surface. 
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Figure 2. -Vortex sheet with bound vortices on y-axis and 
distribution of circulation equal to * A 




Figure 3 - Areas of integration for Infinite span wing. 
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(a) Point P ahead of trailing edge (b) Point P behind trailing edge 

Figure 4- Areas of integration for arbitrary plan form with 
supersonic trailing edge. 


Vo 

V 



Figure 5. - Triangular wing swept behind Mach cone 
showing location of regions A and B. 
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Figure 6 . - Variation of the port of downvrash on x-axis 
Induced by doublets in wake with distance downstream 


in chord lengths , x a - x/c a . Triangular wing. 



Figure 7. - Variation of the part of downwosh on x-axis induced 
by doublets on the plan form with distance downstream in 
chord lengths, x 0 = x/c 0 . Triangular wing. 
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Figure 8 - Sketch of velocity vectors of the air before reaching , 
on, and after leaving supersonic airfoil . 



Figure 9. -Downwash factor induced by doublets in wake ,w w E 0 /w 0 , 
plotted against factor representing distance downstream 
from trailing edge, (x a - !)/ & e . Triangular wing. 




Figure iQ,~ VoriutiOn Gf the tOtui uQWriwQSh on X - uXiS behind O inuuyUiuf winy swept 


behind Mach cone with distance downstream in chord lengths, x- x/c* . 
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exact 

approximate 



Figure //. -Downwash at a large distance behind a wing with 
elliptical span loading in either subsonic or supersonic 
flight compared with first order expansion determined 
from equation (44). 



